The  Astrophysical  Journal,  726:16  (13pp),  2011  January  1 

©  2011.  The  American  Astronomical  Society.  All  rights  reserved.  Printed  in  the  U.S.A. 


doi:10.1088/0004-637X/726/l/16 


CHARACTERISTICS  OF  GAMMA-RAY  LOUD  BLAZARS  IN  THE  VLBA  IMAGING  AND 

POLARIMETRY  SURVEY 

J.  D.  Linford1,  G.  B.  Taylor1,  R.  W.  Romani2,  S.  E.  Healey2,  J.  F.  Helmboldt3,  A.  C.  S.  Readhead4,  R.  Reeves4, 

J.  L.  Richards4,  and  G.  Cotter5 

1  Department  of  Physics  and  Astronomy,  University  of  New  Mexico,  MSC07  4220,  Albuquerque,  NM  87131-0001,  USA 
2  Department  of  Physics,  Stanford  University,  Stanford,  CA  94305,  USA 
3  Naval  Research  Laboratory,  Code  7213,  4555  Overlook  Avenue  SW,  Washington,  DC  20375,  USA 
4  Astronomy  Department,  California  Institute  of  Technology,  Mail  Code  247-17,  1200  East  California  Boulevard,  Pasadena,  CA  91125,  USA 
5  Department  of  Astrophysics,  University  of  Oxford,  Denys  Wilkinson  Building,  Keble  Road,  Oxford  OX1  3RH,  UK 
Received  2010  August  4;  accepted  2010  October  27;  published  2010  December  9 


ABSTRACT 

The  radio  properties  of  blazars  detected  by  the  Large  Area  Telescope  (LAT)  on  board  the  Fermi  Gamma-ray  Space 
Telescope  have  been  observed  as  part  of  the  VLBA  Imaging  and  Polarimetry  Survey.  This  large,  flux-limited  sample 
of  active  galactic  nuclei  (AGNs)  provides  insights  into  the  mechanism  that  produces  strong  y-ray  emission.  At 
lower  flux  levels,  radio  flux  density  does  not  directly  correlate  with  y-ray  flux.  We  find  that  the  LAT-detected 
BL  Lac  objects  tend  to  be  similar  to  the  non-LAT  BL  Lac  objects,  but  that  the  LAT-detected  FSRQs  are  often 
significantly  different  from  the  non-LAT  FSRQs.  The  differences  between  the  y-ray  loud  and  quiet  FSRQs  can  be 
explained  by  Doppler  boosting;  these  objects  appear  to  require  larger  Doppler  factors  than  those  of  the  BL  Lac 
objects.  It  is  possible  that  the  y-ray  loud  FSRQs  are  fundamentally  different  from  the  y-ray  quiet  FSRQs.  Strong 
polarization  at  the  base  of  the  jet  appears  to  be  a  signature  for  y-ray  loud  AGNs. 

Key  words:  galaxies:  active  -  galaxies:  jets  -  galaxies:  nuclei  -  gamma  rays:  galaxies  -  radio  continuum:  galaxies 


1.  INTRODUCTION 

Several  outstanding  questions  about  y-ray  loud  active  galactic 
nuclei  (AGNs)  remain,  including:  “why  are  some  sources  y-ray 
loud  and  others  are  y-ray  quiet?”  The  answer  to  this  is  likely 
related  to  Doppler  boosting  (Jorstad  et  al.  2001;  Taylor  et  al. 
2007;  Kovalev  et  al.  2009;  Lister  et  al.  2009a;  Savolainen  et  al. 
2010),  though  many  details  remain  to  be  filled  in.  A  related 
question  is  “where  are  the  y- rays  produced:  in  or  near  the  core 
or  throughout  the  jets?”  There  are  arguments  for  both  locations 
(Sikora  et  al.  2009;  Tavecchio  et  al.  2010;  Marscher  2006; 
Marscher  et  al.  2008,  2010).  In  either  case  it  still  is  unclear  how 
the  radio  flaring,  or  the  ejection  of  new  components,  is  related 
to  the  gamma-ray  activity,  or  whether  there  are  multiple  classes 
of  y-ray  blazars  (e.g.,  Sikora  et  al.  2001,  2002).  In  addition, 
Stawarz  et  al.  (2008)  argue  that  we  should  detect  y-rays  from 
the  compact  (<kpc)  lobes  of  young  radio  sources. 

The  Large  Area  Telescope  (LAT;  Atwood  et  al.  2009)  on 
board  the  Fermi  Gamma-ray  Space  Telescope  is  a  wide-field 
telescope  covering  the  energy  range  from  about  20  MeV  to 
more  than  300  GeV.  At  high  galactic  latitudes  80%  (106  of  132) 
of  the  y-ray  bright  sources  detected  in  the  LAT  Bright  Source 
List  derived  from  the  first  three  months  of  Fermi  observations 
(Abdo  et  al.  2009a)  are  associated  with  known  AGNs.  This 
trend  continues  in  the  LAT  first-year  catalog  (1FGL;  Abdo  et  al. 
2010a),  where  662  of  1043  high  latitude  sources  are  associated 
with  AGNs  with  high  confidence  (P  >  80%).  The  vast  majority 
of  the  LAT  y-ray  sources  are  blazars,  with  strong,  compact 
radio  emission.  These  blazars  exhibit  flat  radio  spectra,  rapid 
variability,  compact  cores  with  one-sided  parsec-scale  jets,  and 
superluminal  motion  in  the  jets  (Marscher  2006).  A  handful 
of  radio  galaxies  and  Seyferts  have  also  been  identified  in  the 
First  LAT  AGN  Catalog  (1LAC;  Abdo  et  al.  2010b).  These 
radio  galaxies  (e.g.,  3C84,  M87,  Cen  A)  have  faint  parsec-scale 
counterjets  indicating  that  they  are  viewed  at  larger  angles  to 
the  line  of  sight  than  the  blazars.  Also,  there  is  a  much  larger 


fraction  of  unclassified  sources  in  the  one  year  catalog  than  was 
present  in  the  three-month  catalog. 

The  VLBA  Imaging  and  Polarimetry  Survey  (VIPS; 
Helmboldt  et  al.  2007)  was  originally  devised  to  obtain  Very 
Long  Baseline  Array  (VLBA)  measurements  of  AGNs  in  the 
fifth  data  release  of  the  Sloan  Digital  Sky  Survey  (DR5; 
Adelman-McCarthy  et  al.  2006).  It  consists  of  5  GHz  VLBA 
images  and  polarization  data  of  1 127  AGNs  located  between 
declinations  of  15°  and  65°.  The  VIPS  sources  were  drawn  from 
the  Cosmic  Lens  All-Sky  Survey  (CLASS;  Myers  et  al.  2003). 
All  VIPS  sources  were  required  to  have  CLASS  flux  densities 
of  85  mJy  or  more.  The  VLBA  at  5  GHz  gives  us  an  angular  res¬ 
olution  of  just  under  2  mas,  which  allows  us  to  probe  sub-parsec 
scale  features  of  the  objects.  Using  multi-epoch  observations, 
we  will  be  able  to  measure  the  apparent  jet  velocities. 

Various  programs,  such  as  MOJAVE  observing  at  15  GHz 
(Lister  et  al.  2009a;  Homan  et  al.  2009),  the  Boston  University 
program  observing  at  22  and  43  GHz  (Roca-Sogorb  et  al.  2010), 
and  TANAMI  (Ojha  et  al.  2010)  observing  at  8.4  and  22  GHz, 
are  in  place  to  monitor  the  radio  jets  from  the  brightest  blazars 
such  as  3C273,  BL  Lacertae,  etc.  The  detection  of  Centaurus  A, 
3C84  and  M87  by  Fermi  (Abdo  et  al.  2009a,  2009b,  2009c)  and 
the  detection  of  M87  at  TeV  energies  (Aharonian  et  al.  2006) 
suggests  that  other  classes  of  AGNs  (e.g.,  radio  galaxies)  may 
form  a  significant  population  of  y-ray  emitters.  Many  of  the 
current  programs  are  strongly  biased  toward  the  radio  bright 
sample.  VIPS  has  a  flux  limit  roughly  an  order  of  magnitude 
below  the  MOJAVE  survey  and  most  other  samples  that  have 
been  used  to  study  the  properties  of  y-ray  blazars  and  so  allows 
us  to  probe  the  extension  of  the  radio  core/y-ray  properties  down 
to  a  fainter  population.  Many  different  radio-y-ray  correlations 
have  been  suggested  (Taylor  et  al.  2007;  Abdo  et  al.  2009a; 
Lister  et  al.  2009a;  Kovalev  et  al.  2009;  Giroletti  et  al.  2010). 
By  examining  a  larger  sample  (1127  objects)  we  attempt  to 
obtain  more  definitive  insight  into  the  properties  of  y  -ray  loud 
blazars. 
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Figure  1.  CLASS  8.5  GHz  radio  flux  density  (in  mJy)  vs.  redshift  of  LAT- 
detected  sources.  BL  Lac  objects  are  the  top  plot,  FSRQs  are  the  bottom  plot. 


In  Section  2,  we  define  our  sample.  In  Section  3,  we  present 
data  on  the  y-ray  and  radio  properties  of  the  LAT  sources.  In 
Section  4,  we  compare  several  parameters  of  the  LAT  and  non- 
LAT  sources.  Notes  on  some  individual  objects  of  interest  are 
presented  in  Section  5  and  in  Section  6  we  discuss  our  results. 
Some  conclusions  are  given  in  Section  7.  Throughout  this  paper 
we  assume  Ho  —  71  km  s_1  Mpc~*  and  ACDM  cosmology 
(e.g.,  Hinshaw  et  al.  2009). 

2.  SAMPLE  DEFINITION 

We  have  compared  the  11  month  1LAC  to  the  VIPS 
(Helmboldt  et  al.  2007),  and  from  the  intersection  compiled 
a  sample  of  102  LAT  sources  at  \b\  >  10°  associated  with  an 
AGN  with  a  probability  of  more  than  50%  (although  most  of  the 
probabilities  are  greater  than  98%).  All  of  these  objects  were 
observed  with  the  VLBA  from  1996  to  2006  (although  958  of 
the  1 127  sources  were  observed  during  2006)  at  5  GHz.  Each 
VIPS  source  was  observed  for  approximately  500  s,  with  scans 
spread  out  over  1 1  hr.  See  Table  1  for  a  summary  of  our  data.  Of 
the  LAT  objects,  40  are  BL  Lac  objects,  50  are  Flat  Spectrum 
Radio  Quasars  (FSRQs),  and  12  are  either  radio  galaxies  or  un¬ 
classified  objects.  The  optical  classifications  were  adopted  from 
the  Candidate  Gamma-Ray  Blazar  Survey  (CGRaBS;  Healey 
et  al.  2008),  except  that  objects  labeled  in  CGRaBS  as  “AGN” 
we  call  radio  galaxies. 

To  build  a  y-ray  quiet  sample  for  comparison,  we  excluded 
all  LAT  sources  from  the  VIPS  sample,  leaving  1018  objects 
in  our  non-LAT  sample.  Of  these  1018  objects,  24  are  BL  Lac 
objects,  479  are  FSRQs,  and  515  are  either  radio  galaxies  or 
unclassified  objects. 

3.  GAMMA-RAY  FLUX  AND  RADIO  FLUX  DENSITY 

The  radio  flux  densities  in  this  section,  and  those  listed  in 
Table  1,  are  measured  at  8.5  GHz  from  CLASS  (Myers  et  al. 
2003).  For  this  section  we  also  include  seven  MOJAVE-detected 
y- ray  loud  sources  for  which  we  had  CLASS  flux  densities, 
increasing  the  size  of  our  LAT  sample  to  109  sources. 

3.1.  Redshift  Selection  Effect 

Throughout  the  rest  of  this  paper,  we  will  use  the  nonpara- 
metric  Spearman  test  (e.g..  Press  et  al.  1986)  to  look  for  cor¬ 


Figure  2.  LAT  y-ray  flux  (100  MeV-100  GeV)  vs.  redshift.  The  y-ray  fluxes 
are  in  units  of  10  9  photons  cm~2  s'  .  BL  Lac  objects  are  the  top  plot,  FSRQs 
are  the  bottom  plot. 


relations  between  the  LAT-detected  and  non-LAT-detected  ob¬ 
jects.  The  Spearman  correlation  coefficient  (ps)  has  a  range  of 
0  <  | Pi  |  <  1.A  high  value  of  ps  indicates  a  significant  cor¬ 
relation.  This  is  a  powerful  test  for  statistical  correlation,  but  it 
does  not  test  an  actual  physical  correlation.  In  order  to  be  certain 
our  correlations  are  physically  significant,  we  make  sure  that  a 
redshift  selection  effect  is  not  adding  a  bias  to  our  data. 

Imagine  a  population  of  sources  covering  a  wide  range  of 
redshifts  in  which  the  radio  and  y- ray  emission  is  not  physically 
correlated,  but  in  which  there  is  a  correlation  of  radio  flux,  and 
gamma-ray  flux,  with  redshift.  This  is  actually  what  one  would 
naively  expect  because  the  more  distant  objects  will  be  fainter. 
Such  a  population  will  show  a  strong  correlation  of  radio  flux 
with  gamma-ray  flux  via  the  Spearman  test,  but  this  does  not 
indicate  a  significant  physical  correlation  between  these  two 
observables. 

However,  by  investigating  the  relationship  between  radio  flux 
density  and  redshift  for  the  LAT-detected  sources  (see  Figure  1) 
we  can  rule  out  a  correlation  between  the  two.  We  calculated  the 
Spearman  ps  values  for  the  .SY 5-"  relationship  for  the  BL  Lac 
objects  and  FSRQs  separately.  For  the  BL  Lac  objects,  the  ps 
value  is  0.0737,  with  a  79%  probability  that  random  sampling 
would  produce  this  same  ps  value.  For  the  FSRQs,  the  ps  value 
is  0.0450,  with  a  74%  probability  of  getting  the  same  value  by 
random  sampling.  Therefore,  there  is  no  significant  correlation 
between  the  radio  flux  density  and  redshift  for  our  objects. 

For  completeness,  we  also  tested  the  correlation  between 
y-ray  flux  and  redshift.  The  ps  values  are  —0.319  for  BL  Lac 
objects  and  0.104  for  FSRQs,  with  the  probability  of  getting 
the  same  values  from  random  sampling  of  21%  and  45%, 
respectively.  So,  there  is  no  significant  correlation  between 
y-ray  flux  and  redshift,  either.  See  Figure  2  for  a  plot  of  y-ray 
flux  versus  redshift. 

3.2.  Gamma-ray  Flux  versus  Radio  Flux  Density 

In  Figure  3,  we  plot  the  LAT  flux  versus  the  flux  density 
at  8.5  GHz  from  the  CLASS.  The  LAT  fluxes  are  broadband 
fluxes  from  100  MeV  to  100  GeV.  The  BL  Lac  objects  have  a 
large  range  in  their  flux  density  and  y-ray  flux,  but  one  can  see 
a  large  fraction  of  them  have  relatively  low  radio  flux  density 
(5g  5  <  400  mJy)  but  high  y-ray  flux.  The  FSRQ  objects  also 
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Table  1 

LAT  Sources  in  VIPS 


1FGL  Name 

VIPS  Name 

R.A. 

Decl. 

Probability 

Type 

Z 

F35 

A  F35 

Core  Tg 

Opening  Angle 

AP.A. 

5(8.5  GHz) 

1FGL  J0742.2+5443 

J07426+5444 

07:42:39.8 

+54:44:24.679 

i 

FSRQ 

0.723 

2.2 

0.3 

5.26E+10 

142.8 

1FGL  J0746. 6+2548 

J07464+2549 

07:46:25.9 

+25:49:02.146 

0.98 

FSRQ 

2.979 

0.7 

0.2 

5.07E+1 1 

291.7 

1FGL  J0752.8+5353 

J07530+5352 

07:53:01.4 

+53:52:59.636 

1 

BLL 

0.2 

1.2 

0.3 

6.79E+1 1 

1196.9 

1FGL  J0800.5+4407 

J080 11+4401 

08:01:08.3 

+44:01:10.148 

0.92 

0.7 

0.3 

7.36E+10 

230.6 

1FGL  J0806. 2+6148 

J08053+6144 

08:05:18.2 

+61:44:23.704 

0.94 

FSRQ 

3.033 

0.6 

0.2 

1.30E+11 

725.7 

1FGL  J0809.4+3455 

J08096+3455 

08:09:38.9 

+34:55:37.248 

0.99 

BLL 

0.082 

0.7 

0.3 

2.48E+10 

152.1 

1FGL  J0809.5+5219 

J08098+5218 

08:09:49.2 

+52:18:58.252 

0.99 

BLL 

0.138 

2.4 

0.3 

4.83E+09 

154.2 

1FGL  J0815. 0+6434 

J08 146+6431 

08:14:39.2 

+64:31:22.040 

0.99 

BLL 

1.7 

0.3 

2.79E+10 

221.6 

1FGLJ08 16.7+5739 

J08163+5739 

08:16:23.8 

+57:39:09.509 

1 

BLL 

0.7 

0.3 

4.16E+07 

52 

8 

92 

1FGL  J0818.2+4222 

J08 182+4222 

08:18:16.0 

+42:22:45.408 

1 

BLL 

8.7 

0.6 

3.76E+10 

1046.4 

1FGL  J0825.0+5555 

J08247+5552 

08:24:47.2 

+55:52:42.662 

0.99 

FSRQ 

1.417 

0.9 

0.3 

7.86E+10 

15 

-7 

1667.8 

1FGL  J0834.4+4221 

J08338+4224 

08:33:53.9 

+42:24:01.859 

0.98 

FSRQ 

0.249 

0.9 

0.2 

8.66E+09 

560.8 

1FGL  J0856.6+2103 

J08566+2057 

08:56:39.7 

+20:57:43.426 

0.94 

BLL 

1 

0.3 

8.98E+09 

132.7 

1FGL  J0856.6+2103 

J08569+2111 

08:56:57.2 

+21:11:43.640 

0.96 

FSRQ 

2.098 

1 

0.3 

1.70E+10 

332.8 

1FGLJ09 10.7+3332 

J09106+3329 

09:10:37.0 

+33:29:24.418 

1 

BLL 

0.354 

0.8 

0.3 

7.40E+10 

88.6 

1FGL  J0911.0+2247 

J09 107+2248 

09:10:42.1 

+22:48:35.565 

0.98 

FSRQ 

2.661 

2 

0.3 

8.82E+10 

178.8 

1FGL  J0912. 3+4127 

J09121+4126 

09:12:11.6 

+41:26:09.356 

0.98 

FSRQ 

2.563 

0.6 

0.2 

8.06E+09 

166.9 

1FGL  J0915.7+2931 

J09158+2933 

09:15:52.4 

+29:33:23.982 

1 

BLL 

2.1 

0.3 

1.87E+10 

37 

6 

186.1 

1FGL  J0919.6+6216 

J092 16+62 15 

09:21:36.2 

+62:15:52.185 

0.82 

FSRQ 

1.446 

1.1 

0.3 

2.61E+1 1 

1532 

1FGL  J0920.9+4441 

J09209+4441 

09:20:58.5 

+44:41:53.988 

1 

FSRQ 

2.19 

14 

0.7 

8.58E+10 

1368.3 

1FGL  J0923.2+4121 

J09235+4125 

09:23:31.3 

+41:25:27.429 

0.97 

RG 

0.028 

1.2 

0.3 

4.91E+09 

27 

-8 

235.2 

1FGL  J0924.2+2812 

J09238+2815 

09:23:51.5 

+28:15:24.966 

0.97 

FSRQ 

0.744 

1.3 

0.3 

9.87E+1 1 

218.4 

1FGL  J0929. 4+5000 

J09292+5013 

09:29:15.4 

+50:13:35.982 

0.67 

BLL 

0.7 

0.2 

6.67E+10 

747.8 

1FGL  J0934.5+3929 

J0934 1+3926 

09:34:06.7 

+39:26:32.125 

0.93 

BLL 

0.6 

0.2 

1.51E+10 

188.3 

1FGL  J0937.7+5005 

J09372+5008 

09:37:12.3 

+50:08:52.082 

0.97 

FSRQ 

0.276 

0.8 

0.2 

6.26E+10 

372.5 

1FGLJ094 1.2+2722 

J094 18+2728 

09:41:48.1 

+27:28:38.818 

0.97 

FSRQ 

1.306 

0.9 

0 

1.94E+10 

268.4 

1FGL  J0949. 8+1757 

J09496+1752 

09:49:39.8 

+17:52:49.432 

0.98 

FSRQ 

0.693 

0.9 

0 

1.27E+1 1 

153 

1FGL  J0956.9+25 1 3 

J09568+2515 

09:56:49.9 

+25:15:16.047 

0.99 

FSRQ 

0.712 

0.7 

0.2 

1.55E+10 

1905 

1FGL  J0957.7+5523 

J09576+5522 

09:57:38.2 

+55:22:57.740 

1 

FSRQ 

0.896 

10.5 

0.6 

1.20E+08 

60 

-33 

1498.9 

1FGL  J1000.9+2915 

J1001 1+2911 

10:01:10.2 

+29:11:37.543 

0.99 

BLL 

1.4 

0.3 

3.59E+10 

262.1 

1FGL  J1015. 1+4927 

J101 50+4926 

10:15:04.1 

+49:26:00.704 

1 

BLL 

0.2 

8.7 

0.6 

8.16E+09 

16 

-1 

252.2 

1FGL  J1023.6+3937 

J 1023 1+3948 

10:23:11.6 

+39:48:15.378 

0.89 

FSRQ 

1.254 

0.6 

0.2 

9.53E+10 

16 

9 

779.6 

1FGL  J1033.2+4116 

J10330+41 16 

10:33:03.7 

+41:16:06.234 

0.99 

FSRQ 

1.117 

1.1 

0.3 

4.84E+10 

49 

163 

381.4 

1FGL  J1033. 8+6048 

J10338+6051 

10:33:51.4 

+60:51:07.342 

1 

FSRQ 

1.401 

2.2 

0.3 

1.16E+10 

430.4 

1FGL  J1043. 1+2404 

J 1043 1+2408 

10:43:09.0 

+24:08:35.430 

0.99 

BLL 

0.56 

1.5 

0.3 

1.56E+10 

685.4 

1FGL  J1058.6+5628 

J 105 86+5628 

10:58:37.7 

+56:28:11.180 

1 

BLL 

0.143 

5.7 

0.5 

9.49E+09 

189.8 

1FGLJ1 104.4+3812 

J1 1044+3812 

11:04:27.3 

+38:12:31.794 

1 

BLL 

0.03 

26.1 

1 

2.13E+1 1 

15 

-1 

631.6 

1FGLJ1 106.5+2809 

J1 1061+2812 

11:06:07.3 

+28:12:47.045 

0.97 

0.9 

0.2 

5.08E+10 

193.2 

1FGLJ1 112.8+3444 

J1 1126+3446 

11:12:38.8 

+34:46:39.124 

0.99 

FSRQ 

1.949 

0.9 

0.3 

6.95E+10 

160.5 

1FGLJ1 123.9+2339 

J1 1240+2336 

11:24:02.7 

+23:36:45.876 

1 

BLL 

0.9 

0 

1.07E+10 

472.7 

1FGLJ1 141.8+1549 

J1 1421+1547 

11:42:07.7 

+15:47:54.202 

0.99 

1 

0.3 

5.62E+10 

20 

-15 

274 
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1FGL  Name 

VIPS  Name 

R.A. 

Decl. 

Probability 

Type 

Z 

F35 

A  F35 

Core  Tb 

Opening  Angle 

AP.A. 

S  (8.5  GHz) 

1FGLJ  1146.8+4004 

J1 1469+3958 

11:46:58.3 

+39:58:34.307 

0.98 

FSRQ 

1.089 

1 

0.3 

4.80E+1 1 

571.5 

1FGLJ  1150.2+2419 

J1 1503+2417 

11:50:19.2 

+24:17:53.852 

1 

BLL 

0.2 

1 

0.3 

7.60E+10 

17 

7 

745.7 

1FGLJ1 15 1.6+5857 

J1 1514+5859 

11:51:24.7 

+58:59:17.552 

1 

BLL 

0.7 

0.3 

3.08E+09 

35 

10 

102.2 

1FGLJ1 152.1+6027 

J1 1540+6022 

11:54:04.5 

+60:22:20.785 

0.73 

0.6 

0.3 

1.71E+1 1 

151.6 

1FGL  J1202.9+6032 

J 12030+6031 

12:03:03.5 

+60:31:19.129 

1 

RG 

0.065 

1.1 

0.3 

7.22E+09 

21 

12 

191.4 

1FGL  J1209.3+5444 

J 12089+5441 

12:08:54.3 

+54:41:58.190 

0.97 

FSRQ 

1.344 

0.6 

0.2 

1.17E+1 1 

262 

1FGL  J1209.4+4119 

J1 2093+41 19 

12:09:22.8 

+41:19:41.360 

1 

BLL 

0.5 

0.2 

2.99E+10 

485.1 

1FGL  J1209.7+1806 

J12098+1810 

12:09:51.8 

+18:10:06.796 

0.97 

FSRQ 

0.845 

0.8 

0.2 

4.41E+10 

136.3 

1FGL  J1217. 7+3007 

J1 2178+3007 

12:17:52.1 

+30:07:00.625 

1 

BLL 

0.13 

6.7 

0.6 

5.00E+10 

19 

0 

334.9 

1FGLJ  1220.2+3432 

J 1220 1+3431 

12:20:08.3 

+34:31:21.711 

1 

BLL 

0.6 

0.2 

3.05E+10 

314.1 

1FGL  J1221 .5+2814 

J12215+2813 

12:21:31.7 

+28:13:58.497 

1 

BLL 

0.102 

6.9 

0.5 

2.30E+10 

27 

-12 

893.1 

1FGL  J1225. 8+4336 

J1 2248+4335 

12:24:51.5 

+43:35:19.276 

0.94 

0.8 

0 

1.80E+10 

220.8 

1FGL  J1225. 8+4336 

J 12269+4340 

12:26:57.9 

+43:40:58.438 

0.87 

FSRQ 

2.002 

0.8 

0 

1.07E+10 

145.1 

1FGL  J1228.2+4855 

J12288+4858 

12:28:51.8 

+48:58:01.293 

0.91 

FSRQ 

1.722 

1.3 

0.3 

1.89E+12 

17 

-6 

265.5 

1FGL  J1230.4+2520 

J 12302+25 18 

12:30:14.1 

+25:18:07.145 

0.99 

BLL 

0.135 

1.4 

0.3 

1.96E+10 

30 

-9 

326.4 

1FGL  J1248.2+5820 

J1 2483+5820 

12:48:18.8 

+58:20:28.725 

1 

BLL 

4.5 

0.4 

1.54E+10 

43 

15 

324.4 

1FGL  J1253.0+5301 

J12531+5301 

12:53:11.9 

+53:01:11.741 

1 

BLL 

3 

0.4 

3.05E+10 

14 

5 

366 

1FGL  J1258. 3+3227 

J1 2579+3229 

12:57:57.2 

+32:29:29.321 

0.98 

FSRQ 

0.806 

0.8 

0 

2.26E+1 1 

477.1 

1FGLJ1 303.0+2433 

J1 3030+2433 

13:03:03.2 

+24:33:55.684 

1 

BLL 

3.5 

0.4 

7.16E+10 

135.7 

1FGL  J1308.5+3550 

J 13083+3546 

13:08:23.7 

+35:46:37.160 

0.99 

FSRQ 

1.055 

2.2 

0.3 

2.15E+1 1 

60 

-4 

527 

1FGLJ13 12.4+4827 

J13 127+4828 

13:12:43.4 

+48:28:30.928 

0.99 

FSRQ 

0.501 

1.4 

0.3 

2.24E+10 

250.3 

1FGL  J1314.7+2346 

J13147+2348 

13:14:43.8 

+23:48:26.701 

1 

BLL 

1.4 

0.3 

6.62E+09 

158.2 

1FGL  J1317. 8+3425 

J13 176+3425 

13:17:36.5 

+34:25:15.921 

0.98 

FSRQ 

1.05 

0.5 

0.2 

1.54E+10 

541.3 

1FGLJ1321. 1+2214 

J1321 1+2216 

13:21:11.2 

+22:16:12.098 

0.99 

FSRQ 

0.943 

1.1 

0.3 

1.83E+10 

323.6 

1FGLJ133 1.0+5202 

J13307+5202 

13:30:42.6 

+52:02:15.448 

0.99 

RG 

0.688 

0.5 

0 

5.43E+10 

171.9 

1FGL  J1332.9+4728 

J13327+4722 

13:32:45.2 

+47:22:22.653 

0.96 

FSRQ 

0.669 

0.5 

0.2 

1.02E+11 

457.5 

1FGL  J1333.2+5056 

J13338+5057 

13:33:53.8 

+50:57:35.914 

0.93 

1.4 

0.3 

6.24E+09 

107.6 

1FGL  J1345.4+4453 

J13455+4452 

13:45:33.2 

+44:52:59.581 

0.99 

FSRQ 

2.534 

1.6 

0.3 

2.01E+1 1 

197.5 

1FGL  J1351.0+3035 

J13508+3034 

13:50:52.7 

+30:34:53.582 

1 

FSRQ 

0.714 

0.6 

0.2 

3.31E+10 

719.9 

1FGL  J1359. 1+5539 

J13590+5544 

13:59:05.7 

+55:44:29.362 

0.98 

FSRQ 

1.014 

0.9 

0.2 

9.56E+10 

143.8 

1FGL  J1421.0+5421 

J14197+5423 

14:19:46.6 

+54:23:14.783 

0.71 

BLL 

0.152 

0.9 

0.2 

5.81E+10 

20 

4 

2247.3 

1FGL  J1426. 9+2347 

J 14270+2348 

14:27:00.4 

+23:48:00.045 

1 

BLL 

10.2 

0.6 

8.28E+09 

73 

-84 

241.6 

1FGL  J1433.9+4204 

J14340+4203 

14:34:05.7 

+42:03:16.010 

1 

FSRQ 

1.24 

0.7 

0.2 

4.29E+10 

286.8 

1FGL  J1436.9+2314 

J14366+2321 

14:36:41.0 

+23:21:03.297 

0.96 

FSRQ 

1.545 

0.5 

0.2 

1.40E+11 

650.7 

1FGL  J1438.7+371 1 

J14388+3710 

14:38:53.6 

+37:10:35.408 

0.99 

FSRQ 

2.401 

0.7 

0.2 

1.45E+09 

35 

22 

315.4 

1FGL  J145 1.0+5204 

J14509+5201 

14:51:00.0 

+52:01:11.700 

0.99 

BLL 

0.9 

0.3 

2.86E+09 

86.9 

1FGL  J1454.6+5125 

J14544+5124 

14:54:27.1 

+51:24:33.734 

1 

BLL 

1.1 

0.3 

6.16E+09 

96.3 

1FGL  J1503. 3+4759 

J 15037+4759 

15:03:48.0 

+47:59:31.024 

0.96 

BLL 

0.6 

0.2 

4.09E+10 

94.1 

1FGL  J1505. 8+3725 

J15061+3730 

15:06:09.5 

+37:30:51.128 

0.98 

FSRQ 

0.674 

0.7 

0.2 

1.40E+1 1 

18 

4 

777.3 

1FGL  J1516.9+1928 

J15169+1932 

15:16:56.8 

+19:32:13.010 

1 

BLL 

0.6 

0.3 

3.94E+10 

609.3 

1FGL  J1522. 1+3143 

J15221+3144 

15:22:10.0 

+31:44:14.427 

1 

FSRQ 

1.487 

15.9 

0.8 

1.45E+1 1 

493.8 
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1FGL  Name 

VIPS  Name 

R.A. 

Decl. 

Probability 

Type 

Z 

F35 

A  F35 

Core  Tg 

Opening  Angle 

AP.A. 

S  (8.5  GHz) 

1FGL  J1539.7+2747 

J1 5396+2744 

15:39:39.1 

+27:44:38.288 

i 

FSRQ 

2.19 

1 

0.2 

1.29E+1 1 

291.5 

1FGL  J1542.9+6129 

J 15429+6 129 

15:42:56.9 

+61:29:55.358 

i 

BLL 

5.2 

0.4 

5.94E+10 

144 

1FGL  J1558.9+5627 

J1 5588+5625 

15:58:48.3 

+56:25:14.108 

i 

BLL 

0.3 

1.8 

0.3 

1.30E+10 

169.7 

1FGL  J1604.3+5710 

J 16046+57 14 

16:04:37.4 

+57:14:36.668 

0.98 

FSRQ 

0.72 

1.2 

0.3 

1.61E+10 

485.3 

1FGL  J1607. 1+1552 

J16071+1551 

16:07:06.4 

+15:51:34.495 

1 

RG 

0.496 

2.4 

0.4 

2.80E+10 

13 

-1 

223.5 

1FGLJ1616. 1+4637 

J1 6160+4632 

16:16:03.8 

+46:32:25.239 

0.96 

FSRQ 

0.95 

0.6 

0.2 

2.77E+1 1 

125.5 

1FGL  J1635.0+3808 

J16352+3808 

16:35:15.5 

+38:08:04.497 

0.99 

FSRQ 

1.814 

6.8 

0.5 

2.47E+10 

2423.3 

1FGL  J1637.9+4707 

J16377+4717 

16:37:45.1 

+47:17:33.822 

0.91 

FSRQ 

0.74 

1 

0 

6.85E+10 

815.6 

1FGL  J1637. 9+4707 

J 16399+4705 

16:39:56.0 

+47:05:23.575 

0.51 

FSRQ 

0.86 

1 

0 

2.28E+09 

108.7 

1FGL  J1647. 4+4948 

J 16475+4950 

16:47:34.9 

+49:50:00.586 

0.99 

RG 

0.047 

1.1 

0.3 

2.47E+10 

189 

1FGL  J1653.9+3945 

J1 6538+3945 

16:53:51.8 

+39:45:31.490 

0.99 

BLL 

0.034 

8.3 

0.6 

8.10E+10 

17 

12 

1153.2 

1FGL  J1656.9+6017 

J16568+6012 

16:56:48.2 

+60:12:16.455 

0.96 

FSRQ 

0.623 

0.8 

0.2 

2.51E+1 1 

226.5 

1FGLJ  1709.6+4320 

J 17096+43 18 

17:09:41.1 

+43:18:44.547 

1 

FSRQ 

1.027 

1.5 

0.3 

9.36E+10 

40 

10 

197.2 

1FGL  J1724. 0+4002 

J 17240+4004 

17:24:05.4 

+40:04:36.457 

1 

RG 

1.049 

2.9 

0.4 

6.08E+10 

9 

7 

296.1 

1FGL  J1727. 3+4525 

J17274+4530 

17:27:27.6 

+45:30:39.743 

0.97 

FSRQ 

0.714 

1.2 

0.3 

1.95E+11 

1360.3 

1FGL  J1727. 9+5010 

J17283+5013 

17:28:18.6 

+50:13:10.480 

0.98 

BLL 

0.055 

0.9 

0.2 

4.82E+09 

18 

22 

161.8 

1FGL  J1734.4+3859 

J17343+3857 

17:34:20.6 

+38:57:51.446 

1 

FSRQ 

0.976 

6 

0.5 

1.52E+11 

1101.3 

1FGL  J1740.0+5209 

J 17406+52 11 

17:40:37.0 

+52:11:43.413 

0.78 

FSRQ 

1.379 

3.4 

0.4 

5.75E+10 

1358.5 

1FGLJ1742. 1+5947 

J17425+5945 

17:42:32.0 

+59:45:06.729 

0.97 

BLL 

0.6 

0.2 

3.86E+10 

116.2 

1FGL  J1749.0+4323 

J 17490+4321 

17:49:00.4 

+43:21:51.287 

1 

BLL 

2.3 

0.3 

6.17E+10 

21 

-10 

284.9 

Notes.  Column  1:  1FGL  source  name.  Column  2:  VIPS  source  name.  Column  3:  right  ascension  (J2000).  Column  4:  declination  (J2000).  Column  5:  probability  that  the  LAT  source  is  associated  with  the  VIPS  source. 
Column  6:  source  classification.  Column  7:  redshift.  Column  8:  y-ray  flux  in  units  of  10-9  photons  cm-2  s-1  for  1-100  GeV.  Note:  this  is  not  the  full  100  MeV-100  GeV  flux.  Column  9:  error  in  the  y-ray  flux  in  the 
same  units  as  the  flux.  Column  10:  core  brightness  temperature  as  measured  by  automated  program  (Helmboldt  et  al.  2007).  Column  11:  opening  angle  in  degrees.  Column  12:  change  in  jet  position  angle  in  degrees. 
Column  13:  CLASS  flux  density  at  8.5  GHz  in  mJy. 
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Radio  flux  density  at  8.5  GHz  (mJy) 


Figure  3.  LAT  y- ray  flux  (100  MeV-100  GeV)  vs.  CLASS  radio  flux  density 
at  8.5  GHz.  The  y- ray  fluxes  are  in  units  of  1 0  9  photons  cm-2  s_1 .  The  black 
squares  are  BL  Lac  objects,  the  blue  triangles  are  FSRQs,  and  the  red  circles 
are  radio  galaxies  and  unclassified  objects. 


have  a  very  large  range  in  both  their  radio  flux  density  and 
their  y  -ray  flux,  but  the  higher  flux  density  FSRQs  also  tend 
to  have  higher  y-ray  flux.  Overall,  there  is  a  lack  of  strong 
y- ray  emitters  among  the  weak  radio  sources,  but  the  weak 
y -ray  emitters  can  be  bright  or  faint  at  radio  wavelengths. 
The  median  y-ray  flux  is  3.5  x  10~8  photons  cm-2  s'1  for 
BL  Lac  objects,  3.7  x  10~8  photons  cm-2  s-1  for  FSRQs, 
and  3.8  x  10~8  photons  cm-2  s-1  for  the  radio  galaxies  and 
unclassified  objects. 

The  Spearman  ps  values  we  found  were  0.217  for  the  BL 
Lac  objects,  0.318  for  the  FSRQs,  and  —0.413  for  the  radio 
galaxies  and  unclassified  objects.  For  the  BL  Lac  objects  and 
radio  galaxies/unclassified  objects,  the  chances  that  random 
samplings  would  result  in  the  same  ps  values  were  greater  than 
16%,  indicating  a  lack  of  correlation  between  radio  flux  density 
and  y- ray  flux.  For  the  FSRQs,  the  chance  of  a  random  sampling 
having  the  same  ps  was  only  1.7%,  which  indicates  that  there 
is  a  significant  correlation  between  8.5  GHz  flux  density  and 
y-ray  flux  for  the  FSRQs.  However,  when  we  exclude  the 
brightest  10%  of  sources  in  the  radio,  we  find  that  the  remaining 
FSRQs  no  longer  correlate.  The  weaker  90%  of  sources  have 
Pi’s  of  0.135  for  BL  Lac  objects  and  0.159  for  FSRQs,  with  the 
chances  of  random  samples  having  the  same  ps  values  of  18% 
for  BL  Lac  objects  and  33%  for  FSRQs. 

The  lack  of  correlation  for  BL  Lac  objects  is  counter  to  what 
has  been  reported  in  other  studies.  Kovalev  et  al.  (2009)  find  a 
definite  correlation  for  their  MOJAVE  sources  at  15  GHz.  Their 
radio  data  were  taken  closer  to  the  LAT  measurements,  and  the 
y-ray  blazars  are  suspected  to  have  higher  radio  flux  densities 
during  or  shortly  after  y-ray  emission  (Pushkarev  et  al.  2010). 
However,  their  sample  is  only  about  1  /3  the  size  of  ours  and 
it  is  biased  to  sources  with  high  radio  flux  density  (>1.5  Jy). 
Very  strong  flux-flux  correlations  are  found  by  Ghirlanda  et  al. 
(2010)  at  20  GHz.  They  have  a  larger  sample  than  ours  with 
230  LAT  sources,  about  half  of  which  are  FSRQs.  However, 
their  sample  is  also  mostly  comprised  of  bright  sources.  More 
than  half  of  their  sources  have  Siq  >  500  mJy,  compared  to  our 
sample  where  about  1  /3  have  Ss.5  >  500  mJy. 

To  see  if  we  could  reproduce  the  correlation  seen  by  others, 
we  applied  the  Spearman  test  to  the  36  sources  with  8.5  GHz 


flux  densities  above  500  mJy.  This  subsample  is  roughly  the 
brightest  third  of  our  overall  sample  and  contains  1 1  BL  Lac 
objects  and  25  FSRQs.  At  first,  we  did  not  break  the  subsample 
into  optical  types.  The  result  is  a  ps  of  0.390.  The  chance  of  a 
random  sampling  having  the  same  ps  value  as  our  strong  source 
subsample  was  only  1.9%,  indicating  that  this  is  a  significant 
correlation.  When  the  subsample  was  broken  up  by  optical  type, 
the  ps’ s  were  0.214  for  BL  Lac  objects  and  0.422  for  FSRQs. 
The  chances  of  random  samples  having  the  same  ps’ s  were 
69%  and  3.7%,  respectively.  So  we  find  that  FSRQs  exhibit  a 
significant  correlation  between  radio  flux  density  and  y-ray  flux 
but  BL  Lac  objects  do  not,  even  when  limited  to  the  strongest 
radio  sources.  Any  sample  that  contains  mostly  high  radio  flux 
density  objects  will  automatically  be  dominated  by  FSRQs  and 
will  likely  show  a  correlation  between  radio  flux  density  and 
y-ray  flux. 

3.3.  Radio  Flux  Density 

The  median  flux  density  at  8.5  GHz  is  247  mJy  for  the 
LAT  BL  Lac  objects  and  276  mJy  for  the  non-LAT  BL  Lac 
objects.  The  LAT  FSRQs  appear  to  have  significantly  higher 
flux  densities  than  the  non-LAT  FSRQs,  with  a  median  of 
377  mJy  for  LAT  FSRQs  compared  to  a  median  of  238  mJy 
for  non-LAT  FSRQs.  The  LAT  radio  galaxies  and  unclassified 
objects  also  have  higher  8.5  GHz  flux  densities  than  the  non- 
LAT  sources.  The  median  flux  density  for  LAT  radio  galaxies/ 
unclassified  objects  is  207  mJy  and  the  median  flux  density  for 
non-LAT  radio  galaxies/unclassified  objects  is  119  mJy.  The 
Kolomogorov-Smirnov  (K-S)  test  probability  that  the  LAT  and 
non-LAT  BL  Lac  objects  belong  to  the  same  parent  population  is 
about  81%.  The  K-S  test  probability  that  the  LAT  and  non-LAT 
FSRQs  are  related  is  only  0.3%.  The  radio  galaxy /unclassified 
objects  K-S  test  probability  is  0.01%. 

4.  MORPHOLOGICAL  COMPARISON 
4.1.  Source  Classes 

VIPS  sources  are  classified  as  point  sources  (PSs),  short  jets 
(SJETs),  long  jets  (LJETs),  complex  (CPLX),  or  compact  sym¬ 
metric  object  candidates  (CSO)  by  the  automatic  classification 
scheme  of  Helmboldt  et  al.  (2007).  See  Table  2  for  the  source 
classifications  for  both  LAT  and  non-LAT  sources  in  VIPS. 
From  the  classifications,  it  appears  that  SJETs  and  PS  BL  Lac 
objects  are  less  likely  to  produce  y-rays.  About  75%  of  the  LAT 
BL  Lac  objects  are  classified  as  LJET,  compared  to  only  45% 
for  the  non-LAT  BL  Lac  objects.  The  major  difference  between 
LAT  and  non-LAT  FSRQs  is  the  lack  of  CSOs  in  the  LAT  sam¬ 
ple.  There  is  also  a  lack  of  CSOs  among  the  LAT  radio  galaxies 
and  objects  without  an  optical  classification. 

4.2.  Polarization 

Polarization  is  detected  from  the  cores  of  nearly  half  (49  of 
102)  of  the  LAT  sources  in  our  sample.  We  define  the  core  as  the 
bright  compact  component  at  the  base  of  the  jet.  We  detected 
polarization  in  the  core  of  17  of  the  40  BL  Lac  objects,  24  of 
the  50  FSRQs,  and  8  of  the  12  radio  galaxies  and  unclassified 
objects.  For  the  non-LAT  sources,  core  polarization  is  detected 
in  only  270  of  the  total  1018  sources.  There  are  10  of  24  non- 
LAT  BL  Lac  objects,  158  of  479  FSRQs,  and  102  of  515  radio 
galaxies  or  unclassified  objects  detected.  The  LAT  FSRQs  and 
radio  galaxies/unclassified  objects  tend  to  have  polarization  in 
their  cores  more  than  their  non-LAT  counterparts.  In  contrast. 
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Table  2 

Source  Classifications 


Linford  et  al. 


LAT/Non-LAT 

Optical  Type 

LJET 

SJET 

PS 

CPLX 

CSO 

Unidentified 

LAT-detected 

BL  Lac  objects 

30  (75%) 

4  (10%) 

4  (10%) 

2  (5%) 

FSRQs 

22  (44%) 

1 1  (22%) 

16(32%) 

1  (2%) 

RG/Other 

8  (67%) 

1  (8%) 

3  (25%) 

Non-LAT-detected 

BL  Lac  objects 

1 1  (46%) 

7  (29%) 

6  (25%) 

FSRQs 

188(39%) 

121  (25%) 

136  (28%) 

2(~1%) 

30  (6%) 

2  (~1%) 

RG/Other 

214  (42%) 

98  (19%) 

111  (21%) 

11  (2%) 

71  (14%) 

10  (2%) 

Note.  LJET,  long  jet;  SJET,  short  jet;  PS,  point  source;  CPLX,  complex;  CSO,  compact  symmetric  object  candidate. 
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Figure  4.  Distributions  of  fractional  polarization  in  the  cores  of  LAT-detected 
FSRQs  and  radio  galaxies/unclassified  objects  (top)  and  non-LAT  FSRQs  and 
radio  galaxies/unclassified  objects  (bottom). 

the  LAT  BL  Lac  objects  do  not  appear  to  be  any  more  likely  to 
have  polarized  cores  than  the  non-LAT  BL  Lac  objects.  The 
distributions  of  fractional  polarization  in  the  cores  for  LAT 
and  non-LAT  FSRQs,  and  radio  galaxies/unclassified  objects 
are  shown  in  Figure  4.  The  K-S  test  probabilities  that  these 
distributions  were  drawn  from  the  same  parent  populations  are 
84.6%  for  the  BL  Lac  objects,  39.4%  for  the  FSRQs,  and  76.9% 
for  the  radio  galaxies/unclassified  objects. 

Our  polarization  results  are  different  than  those  reported  by 
Flovatta  et  al.  (2010)  with  their  MOJAVE  data  at  15  GHz.  They 
found  that  the  median  fractional  polarization  for  the  cores  of 
LAT-detected  sources  was  higher  than  for  non-LAT  sources. 
Their  median  values  for  their  2008-2009  observations  are  2.5% 
for  LAT-detected  sources  and  1.86%  for  non-LAT  sources.  Our 
overall  median  values  are  3.5%  for  LAT  sources  and  4.4%  for 
non-LAT  sources.  However,  our  result  is  not  of  high  significance 
as  the  standard  deviations  are  2.5%  for  LAT  sources  and  3.4%  for 
non-LAT  sources.  It  should  also  be  noted  that  our  VIPS  data  are 
not  contemporaneous  with  LAT  observations.  It  is  possible  that 
polarization  increases  with  y-ray  emission.  In  fact,  the  MOJAVE 
data  indicate  that  the  polarization  levels  were  stronger  during 
the  LAT  detections  than  in  previous  years  (2002-2008),  but  only 
by  a  factor  of  1.2. 

4.3.  Core  Brightness  Temperature 

We  have  obtained  the  brightness  temperatures  from  automatic 
modelfits  to  source  components  as  described  in  Helmboldt 
et  al.  (2007).  For  sources  with  two  or  fewer  components  we 


iog(TB) 

Figure  5.  Distributions  of  core  brightness  temperatures  for  LAT-detected  FSRQs 
and  radio  galaxies/unclassified  objects  (top)  compared  to  the  distributions  of 
core  brightness  temperature  for  non-LAT  FSRQs  and  radio  galaxies/unclassified 
objects  (bottom). 

have  further  refined  our  modelfitting  procedure  by  fitting  to  the 
visibility  data  directly  (Taylor  et  al.  2007).  Only  image  plane 
modelfitting  was  carried  out  in  Helmboldt  et  al.  (2007)  due  to 
the  tendency  for  the  automatic  visibility  modelfitting  to  go  awry 
for  complicated  sources. 

There  is  strong  evidence  that  the  brightness  temperature  for 
the  cores  of  the  LAT-detected  sources  may  be  higher  than 
the  y-ray  quiet  population  (see  the  FSRQ  and  radio  galaxy/ 
unclassified  distributions  in  Figure  5).  The  BL  Lac  objects 
are  the  exception  to  this  trend.  The  median  core  brightness 
temperature  for  BL  Lac  objects  detected  by  LAT  is  about 
2.6  x  101()  K,  while  the  median  for  remaining  non-LAT  BL  Lac 
objects  is  2.7  x  1010  K.  The  K-S  test  probability  for  these  two 
distributions  coming  from  the  same  parent  population  is  93.6%. 
The  median  core  brightness  temperature  for  FSRQ  objects  in  the 
LAT  sample  is  8.2  x  1010  K,  while  the  non-LAT  FSRQ  median 
is  only  about  2.5  x  1010  K.  The  K-S  test  probability  that  the 
two  FSRQ  distributions  come  from  the  same  parent  population 
is  very  low  (0.001%).  The  median  core  brightness  temperature 
for  LAT  radio  galaxies  and  unclassified  objects  is  3.9  x  1010  K, 
and  for  non-LAT  objects  the  median  is  1.0  x  1010  K.  The  K-S 
test  probability  that  these  distributions  are  drawn  from  the  same 
parent  sample  is  only  4.4%. 

Interestingly,  y-ray  flux  does  not  appear  to  have  any  relation 
to  core  brightness  temperature.  We  applied  the  Spearman  test 
to  the  y-ray  versus  core  brightness  temperature  distributions 
and  found  correlation  coefficients  of  0.093  for  BL  Lac  objects. 
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Figure  6.  Distribution  of  opening  angles  for  LAT-detected  BL  Lac  objects, 
FSRQs,  and  radio  galaxies/unclassified  objects  (top)  compared  to  the  distribu¬ 
tion  of  opening  angles  for  non-LAT  objects  (bottom). 


0.111  for  FSRQs,  and  —0. 1 8 1  for  radio  galaxies  and  unclassified 
objects. 


4.4.  Jet  Opening  Angle 

We  have  measured  a  mean  opening  half-angle  by  the  fol¬ 
lowing  procedure  (Taylor  et  al.  2007):  We  measure  the  sepa¬ 
ration  of  each  jet  component  from  the  core  along  the  jet  axis 
(taken  to  be  a  linear  fit  to  the  component  positions)  and  the 
distance  of  each  component  from  the  jet  axis,  i.e.,  x'  and  y'  po¬ 
sitions  in  a  rotated  coordinate  system  with  the  jet  axis  along  the 
x'-axis.  For  each  component,  we  measure  its  extent  from  its  cen¬ 
ter  along  a  line  perpendicular  to  the  jet  axis  using  the  parameters 
of  its  elliptical  fit,  and  then  deconvolve  this  using  the  extent  of 
the  Gaussian  restoring  beam  along  the  same  line.  The  opening 
half-angle  measured  from  each  component  is  then  taken  to  be 


et  al.  (2007)  reported  that  EGRET  sources  tended  to  have  greater 
core  fractions  than  sources  not  detected  by  EGRET.  With  the 
LAT  sources,  we  find  that  this  is  not  the  case.  The  median  core 
fractions  are  0.868  for  LAT  BL  Lac  objects  and  0.896  for  non- 
LAT  BL  Lac  objects,  0.939  for  LAT  FSRQs  and  0.929  for  non- 
LAT  FSRQs,  and  0.889  for  LAT  radio  galaxies/unclassified 
objects  and  0.886  for  non-LAT  radio  galaxies/unclassified 
objects.  The  K-S  test  probabilities  that  the  distributions  were 
drawn  from  the  same  parent  distribution  are  26.4%  for  BL 
Lac  objects,  16.1%  for  FSRQs,  and  23.5%  for  radio  galaxies/ 
unclassified  objects,  respectively. 

4.6.  Jet  Bending,  Length,  and  Brightness  Temperature 

We  can  also  compare  the  amount  of  bending  in  the  jet  on  the 
parsec  scale  for  the  LAT  sources  compared  to  non-LAT  sources. 
Unfortunately,  as  with  the  jet  opening  angle  measurements, 
we  are  limited  to  only  30  LAT  sources  and  the  results  are 
inconclusive. 

We  compared  the  distribution  of  jet  length  in  LAT  sources 
with  non-LAT  sources.  There  does  not  appear  to  be  an  apprecia¬ 
ble  difference  between  the  two  populations  of  LSRQ  an  radio 
galaxies/unclassified  objects.  The  K-S  test  results  are  63.3%  for 
the  FSRQ  objects,  and  64.6%  for  the  radio  galaxy  and  unclassi¬ 
fied  objects.  Interestingly,  the  BLLac  objects  in  the  LAT  sample 
appear  to  have  longer  jets  than  the  non-LAT  BL  Lac  objects. 
The  distributions  for  the  BL  Lac  objects  are  shown  in  Figure  8. 
The  K-S  test  probability  that  the  two  are  drawn  from  the  same 
distribution  is  2.1%.  The  FSRQ  and  AGNs/other  objects  have 
essentially  the  same  distributions. 

Unlike  the  core  brightness  temperatures,  the  jet  brightness 
temperatures  (formally,  the  brightness  temperature  of  the  bright¬ 
est  jet  component)  for  the  LAT  and  non-LAT  distributions  look 
fairly  similar.  The  K-S  test  probabilities  for  the  BL  Lac  object, 
LSRQ,  and  radio  galaxy /other  pairs  of  distributions  are  34.2%, 
92.1%,  and  83.5%,  respectively. 

5.  NOTES  ON  INDIVIDUAL  SOURCES 


j/r  =  arctan[(|/|  +dr)/ |x'|], 

where  dr  is  the  deconvolved  Gaussian  size  perpendicular  to 
the  jet  axis.  After  measuring  this  for  each  jet  component,  we 
average  them  to  get  a  single  value.  This  is  only  done  for  sources 
with  more  than  two  total  components,  (i.e.,  at  least  two  jet 
components).  Unfortunately,  we  can  only  make  quantitative 
estimates  for  30  of  the  LAT  sources;  16  BL  Lac,  9  LSRQ,  and  5 
radio  galaxy  or  unclassified  objects.  The  distributions  of  opening 
angles  are  shown  in  Figure  6.  Although  this  is  based  on  very 
small  statistics  there  is  evidence  that  the  LAT  source  jets  have 
unusually  large  opening  angles.  This  is  especially  true  for  the 
FSRQ  detected  by  LAT.  The  MOJAVE  (Pushkarev  et  al.  2009) 
and  TANAMI  (Ojha  et  al.  2010)  groups  both  reported  larger 
opening  angles  for  y-ray  loud  blazars.  Taylor  et  al.  (2007)  also 
found  evidence  for  larger  opening  angles  in  EGRET  sources. 

We  found  10  LAT-detected  objects  with  opening  angles 
greater  than  30°.  The  large  opening  angle  objects  are  BL  Lac 
objects  and  LSRQs,  five  of  each.  See  Ligure  7  for  contour  maps 
of  these  10  sources. 

4.5.  Core  Fraction 

We  define  the  core  fraction  as  the  ratio  of  flux  density  in  the 
core  compared  to  the  total  flux  density  at  5  GHz.  A  high  core 
fraction  can  be  an  indication  of  a  high  Doppler  factor.  Taylor 


5.1.  Low  Core  Brightness  Temperature 

J09576+5522.  With  a  core  brightness  temperature  of  only 
9.7  x  108  K,  this  FSRQ  has  a  y- ray  flux  of  (106.5  ±  7.4)  x 
1  (V9  photons  cm'2  s'1,  nearly  three  times  the  median  for 
FSRQs.  Also  known  as  [HB89]  0954+556,  it  has  a  redshift 
of  0.896  (Abazajian  et  al.  2004)  and  a  8.5  GHz  flux  density 
is  1.5  Jy.  This  FSRQ  has  been  studied  in  the  UV  (Pian  et  al. 
2005),  radio.  X-ray  (Tavecchio  et  al.  2007),  and  was  detected  in 
the  y-rays  by  EGRET  (Hartman  et  al.  1999).  Tavecchio  et  al. 
(2007)  noted  that  X-ray  flux  was  only  found  at  the  terminal 
portion  of  the  jet.  The  automated  program  initially  reported  the 
core  brightness  temperature  for  this  object  to  be  1 .2  x  108  K.  We 
suspected  program  did  not  generate  an  accurate  model  for  this 
source.  When  we  did  the  model  fitting  by  hand  in  DILMAP,  we 
calculated  a  core  brightness  temperature  of  9.7  x  108  K.  This  is 
the  lowest  core  temperature  of  all  objects  in  our  LAT  sample. 
This  object  is  featured  in  Ligure  7  because  its  jet  opening  angle 
is  about  60° . 

J08163+5739.  The  automated  program  found  that  this  BL  Lac 
object  had  a  core  brightness  temperature  of  only  4.16  x  107  K. 
Again,  when  we  did  the  model  fitting  by  hand  we  calculated 
a  core  brightness  temperature  of  1.6  x  109  K.  Even  with  this 
higher  core  brightness  temperature,  this  object  is  still  the  coolest 
BL  Lac  object  in  our  LAT  sample.  This  object  has  a  y-ray  flux  of 
(15.2  ±  0.3)  x  10~9  photons  cm'2  s'1.  It  is  a  compact  and  weak 
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Figure  7.  Contour  maps  of  the  10  LAT  sources  with  opening  angles  larger  than  30°.  All  data  were  taken  with  the  VLBA  at  5  GHz.  Bottom  contours  are  at  0.5  mJy 
beam-1 ,  except  for  the  first  object  which  has  a  bottom  contour  level  is  1  mJy  beam-1 .  The  restoring  beam  is  shown  as  an  ellipse  in  the  lower  left  corner  of  each  panel. 
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Figure  8.  Distribution  of  jet  lengths  for  LAT-detected  BL  Lac  object  (top) 
compared  to  the  distribution  of  jet  lengths  fornon-LAT  BLLac  objects  (bottom). 

radio  source  with  a  radio  flux  density  of  92  mJy  at  8.5  GHz. 
There  is  no  measured  redshift  for  this  object.  It  is  also  known  as 
SBS  0812+578  and  1RXS  J08 1624.6+573910.  It  has  a  small  jet 
extending  to  the  east.  This  object  is  featured  in  Figure  7  because 
its  jet  opening  angle  is  approximately  52°. 

5.2.  High  Redshift 

J07464+2549.  With  a  redshift  of  2.979  (Abazajian  et  al. 
2004),  this  FSRQ  has  a  y-ray  flux  (42.3  ±  7.9)  x  10~9  pho¬ 
tons  cm-2  s_1 .  Its  core  brightness  temperature  is  5.07  x  101 1  K, 
which  is  well  above  the  median  for  the  LAT-detected  FSRQs. 
It  is  also  known  as  B2  0743+25  and  RX  J0746.4+2549.  It  is  a 
relatively  compact  and  bright  radio  source  with  a  flux  density 
of  731  mJy  at  8.5  GHz  (see  Figure  9).  It  has  a  very  small  jet 
extending  to  the  north.  In  2005,  the  Burst  Alert  Telescope  on 
board  the  Swift  satellite  detected  a  flare  of  15-195  keV  emission 
from  this  object  (Sambruna  et  al.  2006). 

J08053+6144.  With  a  redshift  of  3.033  (Sowards-Emmerd 
et  al.  2005),  this  FSRQ  has  the  highest  redshift  among  the  LAT- 
detected  objects.  It  has  a  y-ray  flux  36.0  ±  8.9  x  10'9  pho¬ 
tons  cm-2  s_1.  Its  core  brightness  temperature  is  1 .30  x  101 1  K, 
which  is  slightly  higher  than  the  median  for  LAT-detected  FS¬ 
RQs.  This  object  is  fairly  bright  in  the  radio  with  a  flux  density  of 
725  mJy  at  8.5  GHz.  It  has  a  short  jet  extending  to  the  southeast 
(see  Figure  9). 

5.3.  Radio  Galaxies 

Our  sample  contains  six  sources  identified  as  radio  galaxies, 
although  in  some  cases  there  is  controversy  in  this  identification 
as  described  below. 

J09235+4125.  This  object  is  also  known  as  B3  0920+416 
and  GB6  J0923+4125.  It  has  a  y- ray  flux  of  (40.6  ±  10.0)  x 
10'9  photons  cm-2  s'1  and  a  redshift  of  0.028  (Falco  et  al. 
1998).  Its  peak  radio  flux  density  is  120  mJy  at  8.5  GHz.  It  has  a 
short  jet  extending  to  the  east  (see  Figure  9).  The  NASA/IPAC 
Extragalactic  Database  (NED)  classifies  this  object  as  an  FSRQ. 
However,  this  is  probably  an  error.  Falco  et  al.  (1998)  call  it  a 
“late-type  galaxy.” 

J12030+6031.  This  object  is  also  known  as  SBS  1200+608 
and  RX  J1203.0+6031.  It  has  a  y- ray  flux  of  (44.8  ±  0.3)  x 
1  (L9  photons  cm-2  s'1  and  a  redshift  of  0.065  (Falco  et  al. 


1998).  Its  peak  radio  flux  density  is  105  mJy  at  8.5  GHz.  It  has 
a  short  jet  extending  to  the  south  (see  Figure  9).  Veron-Cetty 
&  Veron  (2006)  classify  this  object  as  a  low  ionization  nuclear 
emission  region  (LINER)  galaxy,  sometimes  referred  to  as  a 
Seyfert  3. 

J13307+5202.  This  radio  galaxy  is  also  known  as  CGRaBS 
J1330+5202.  It  has  an  upper  limit  to  its  y-ray  flux  of  41.5  x 
10'9  photons  cm-2  s'1  and  a  redshift  of  0.688  (Healey  et  al. 
2008).  It  has  a  peak  radio  flux  density  of  127  mJy  at  8.5  GHz 
and  a  jet  extending  to  the  southwest  (see  Figure  9).  Sowards- 
Emmerd  et  al.  (2005)  list  it  as  an  unknown  optical  type. 

J16071+1551.  This  object  is  also  known  as  [HB89] 
1604+159.  It  is  an  EGRET-detected  source  known  as  3EG 
J1605+1553.  It  has  aLAT  y-ray  flux  of  (33.9±7.7)  x  10“9  pho¬ 
tons  cm-2  s“ 1  and  a  redshift  of  0.496  (Adelman-McCarthy  et  al. 
2008).  It  has  a  peak  radio  flux  density  of  214  mJy  at  8.5  GHz 
and  a  jet  extending  to  the  east  (see  Figure  9).  Veron-Cetty  & 
Veron  (2006)  classify  this  object  as  a  BL  Lac  object. 

J16475+4950.  This  object  is  also  known  as  SBS  1646+499 
and  1RXS  J 164735. 4+495001.  It  has  a  y-ray  flux  of  33.3  ± 
8.5  x  10~9  photons  cm-2  s-1  and  a  redshift  of  0.047  (Falco 
et  al.  1998).  It  has  a  peak  radio  flux  density  of  131  mJy  at 
8.5  GHz.  It  has  a  very  short  jet  extending  to  the  southeast  (see 
Figure  9).  The  spectrum  in  Marc  ha  et  al.  (1996)  indicates  that 
this  object  is  a  Seyfert  1  blazar. 

J17240+4004.  This  object  is  also  known  as  B2  1722+40 
and  RX  J1724.0+4004.  It  has  a  y-ray  flux  of  (47.1  ±  9.0)  x 
1(L9  photons  cm'2  s'1  and  a  redshift  of  1.049  (Vermeulen 
et  al.  1996).  This  is  the  brightest  radio  galaxy  in  our  sample 
with  a  peak  flux  density  of  323  mJy  at  8.5  GHz.  It  has  a  long 
jet  extending  to  the  northwest  (see  Figure  9).  Veron-Cetty  & 
Veron  (2006)  tentatively  classify  this  object  as  a  BL  Lac  object. 
However,  the  spectrum  in  Vermeulen  et  al.  (1996)  does  not 
support  this  classification.  The  spectrum  has  only  narrow  lines, 
so  it  cannot  be  an  FSRQ,  and  the  equivalent  widths  are  greater 
than  5  A,  so  it  cannot  be  a  BL  Lac  object. 

6.  DISCUSSION 

One  of  the  primary  explanations  for  some  AGNs  being  y-ray 
loud  is  Doppler  Boosting.  The  orientation  and  material  speed  of 
the  AGN  jets  have  a  large  effect  on  their  observed  fluxes  due  to 
the  Doppler  effect.  The  kinematic  Doppler  factor  for  a  moving 
source  is  defined  as 

S  =  [F(l  —  /3  cos6»)]_I, 

where  fi  is  the  bulk  velocity  v/c,  T  is  the  Lorentz  factor, 
T  =  (1  —  p2)~1/2,  and  6  is  the  angle  between  the  velocity 
vector  and  the  line  of  sight.  Emissions  from  an  approaching 
source  are  blueshifted  by 

v  —  <5v' 

where  v’  is  the  frequency  in  the  source’s  rest  frame.  The  intensity 
enhancement,  commonly  known  as  the  “Doppler  boosting,”  is 
given  in  Rybicki  &  Lightman  (1979)  as 

4,(10  =  S3/>'). 

Because  LAT  measures  the  flux  across  a  large  bandwidth 
(1-100  GeV),  we  have  to  integrate  over  frequency: 

F  =  J  Ivdv  =  J  S3I'v,Sdv'  =  84Ff. 
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peak=7.298E-01  Jy 


peak=7.249E— 01  Jy 


peak=1 .1 24E-01  Jy 


Figure  9.  Contour  maps  of  the  objects  in  Sections  5.2  and  5.3.  All  data  were  taken  with  the  VLBA  at  5  GHz.  Bottom  contours  are  at  0.5  mJy  beam  1 .  The  restoring 
beam  is  shown  as  an  ellipse  in  the  lower  left  corner  of  each  panel. 


Therefore,  a  small  difference  in  material  velocity  or  orientation 
angle  leads  to  a  large  difference  in  intensity.  See  Urry  & 
Padovani  (1995)  for  a  more  thorough  discussion  of  Doppler 
enhancement. 

6.1.  BL  Lac  Objects 

Almost  two-thirds  of  the  BL  Lac  objects  in  VIPS  are  detected 
by  LAT,  and  the  only  significant  difference  between  the  LAT  and 
non-LAT  populations  is  in  their  jet  lengths.  We  have  shown  that 
the  LAT  BL  Lac  objects  have  longer  jets  than  the  non-LAT  BL 
Lac  objects  (see  Section  4.6).  There  are  no  differences  in  the 
fraction  of  polarized  BL  Lac  objects  or  the  distributions  of  the 
polarization.  The  LAT  and  non-LAT  BL  Lac  objects  have  nearly 
identical  core  brightness  temperature  distributions.  Why  are  the 


non-LAT  BL  Lac  objects  quiet  in  the  /-ray?  It  is  possible  that 
the  non-LAT  BL  Lac  objects  are  fundamentally  different  from 
the  LAT  BL  Lac  objects,  but  it  must  be  in  some  subtle  way.  It 
seems  more  likely  that  all  the  BL  Lac  objects  are  emitting  in  the 
/-ray,  but  some  are  emitting  below  the  threshold  for  the  LAT  to 
detect. 

Initially,  we  suspected  that  the  non-LAT  BL  Lac  objects  might 
have  been  too  far  away  to  detect  their  /-rays.  The  more  distant 
BL  Lac  objects  definitely  have  lower  /-ray  flux,  but  it  is  clear 
that  the  non-LAT  BL  Lac  objects  are  not  at  higher  redshifts 
than  the  LAT  BL  Lac  objects.  While  we  only  have  redshift 
measurements  for  16  of  the  LAT  BL  Lac  objects  and  15  of  the 
non-LAT  BL  Lac  objects,  that  is  still  enough  to  rule  out  redshift 
as  the  cause  of  LAT  non-detection.  The  maximum  measured 
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redshift  for  the  LAT  BL  Lac  objects  is  0.56.  Nearly  half  (11  /24) 
of  the  non-LAT  BL  Lac  objects  have  redshifts  lower  than  0.56. 
Also,  there  is  no  strong  correlation  between  redshift  and  y-ray 
flux  for  the  BL  Lac  objects  (see  Section  3.1). 

The  prevalence  of  LJET  BL  Lac  objects  is  an  interesting 
clue  to  this  puzzle.  The  K-S  test  result  of  2.1%  for  the  jet 
length  distributions  (see  Section  4.6)  is  only  of  modest  statistical 
significance,  but  the  fact  that  a  larger  percentage  of  the  LAT  BL 
Lac  objects  are  classified  as  LJET  than  the  non-LAT  BL  Lac 
objects  (see  Table  2)  supports  it.  One  would  expect  the  jets  to 
appear  shorter  if  they  are  pointed  more  directly  at  us.  However, 
if  the  y-ray  loud  BL  Lac  objects  have  higher  Doppler  factors 
than  their  y- ray  quiet  counterparts,  that  would  indicate  that  the 
LAT  BL  Lac  objects  have  more  powerful  jets  than  the  non-LAT 
BL  Lac  objects.  A  more  powerful  jet  will  appear  longer  even 
when  viewed  at  a  small  angle. 

Because  we  find  that  the  LAT  BL  Lac  objects  have  longer  jets 
than  non-LAT  BL  Lac  objects,  the  orientation  angle  for  LAT  BL 
Lac  objects  cannot  be  smaller  than  non-LAT  BL  Lac  objects.  The 
non-LAT  BL  Lac  objects  could  have  slightly  lower  jet  speeds. 
This  can  explain  why  the  LAT  BL  Lac  objects  have  longer  jets 
and  why  they  are  y- ray  loud.  Unfortunately,  this  theory  is  not 
supported  by  recent  work  by  Lister  et  al.  (2009b),  where  they 
find  that  two  non-LAT  BL  Lac  objects  have  higher  maximum 
jet  material  velocities  than  LAT  BL  Lac  objects.  However,  they 
only  had  10  BL  Lac  objects  in  their  LAT-detected  sample,  so 
this  is  still  an  open  question. 

BL  Lac  objects  are  known  to  have  high  variability  in  their 
flux  densities.  The  non-LAT  BL  Lac  objects  could  currently  be 
in  low  y- ray  luminosity  states,  but  may  enter  states  of  enhanced 
y-ray  activity  at  a  later  time.  It  is  possible  that  the  jet  material 
speed  is  related  to  this  flux  variability.  Based  on  the  number  of 
BL  Lac  objects  in  our  sample,  the  duty  cycle  would  need  to  be 
about  2/3  over  the  course  of  one  year. 

6.2.  FSRQs 

The  LAT  FSRQs  seem  to  be  substantially  different  from  the 
non-LAT  FSRQs.  Only  about  9%  (50/529)  of  the  FSRQs  in 
VIPS  are  y-ray  loud.  This  indicates  that  the  LAT  FSRQs  are 
special  in  some  way.  Recall  from  Section  3.2  that  y-ray  flux 
correlates  with  radio  flux  density  for  FSRQs,  especially  for 
high  flux  density  sources.  Also,  in  general  the  LAT  FSRQs 
have  significantly  higher  core  brightness  temperatures,  with  the 
notable  exception  of  J09576+5522,  and  they  appear  to  have 
larger  opening  angles.  So,  brighter  FSRQs  are  more  likely  to 
produce  y-rays.  We  showed  in  Section  3.1  that  y-ray  flux  does 
not  correlate  with  redshift,  so  we  do  not  expect  the  LAT  FSRQs 
to  have  lower  redshift  than  their  non-LAT  counterparts.  We 
explored  the  possibility  that  LAT  FSRQs  were  at  low  redshift, 
but  found  that  they  cover  a  range  similar  to  the  non-LAT  FSRQs. 

The  y-ray  flux  from  FSRQs  can  be  explained  with  Doppler 
boosting,  but  requires  a  more  dramatic  difference  than  with  the 
BL  Lac  objects.  Lister  et  al.  (2009b)  reported  that  LAT-detected 
FSRQs  had  higher  median  jet  speeds  than  LAT-detected  BL 
Lac  objects  by  more  than  a  factor  of  two.  It  is  also  possible  that 
there  are  two  subclasses  for  FSRQs:  y-ray  loud  and  y-ray  quiet. 
The  difference  in  y-ray  detection  rates  for  FSRQs  could  also 
be  due  to  variations  in  the  y-ray  emission,  however,  the  duty 
cycle  required  (about  10%)  is  much  more  extreme  than  for  the 
BL  Lac  objects. 

A  larger  fraction  of  the  LAT  FSRQs  have  polarization  in 
their  cores  compared  to  the  non-LAT  FSRQs.  We  detected 
polarization  in  the  cores  of  48%  of  the  LAT  FSRQs  compared 


to  33%  of  the  non-LAT  FSRQs,  a  difference  that  is  significant  at 
about  the  la  level.  This  points  to  an  intrinsic  difference  in  y-ray 
quiet  versus  y-ray  loud  FSRQs.  The  higher  polarization  could 
indicate  more  organized  magnetic  fields,  possibly  pointing  to 
a  greater  spin  of  the  supermassive  black  hole  (Meier  2001). 
A  higher  black  hole  spin  would  also  lead  to  higher  material 
velocities  in  the  jets  (Meier  1999). 

6.3.  Radio  Galaxies 

There  is  a  definite  lack  of  CSOs  among  the  LAT-detected 
objects.  This  is  not  surprising  from  a  Doppler  boosting  point 
of  view.  For  an  object  to  appear  symmetric,  we  need  to  see 
it  close  to  “edge  on”  so  we  can  see  both  jets.  With  the  jets 
pointed  nearly  perpendicular  to  our  line  of  sight,  we  would  not 
expect  to  see  much  Doppler  boosting  in  the  radiation.  However, 
Stawarz  et  al.  (2008)  suggested  that  CSOs  should  account  for  a 
“relatively  numerous  class  of  extragalactic  sources”  for  Fermi 
due  to  ultrarelativistic  electrons  in  the  lobes  of  these  compact 
objects.  We  have  yet  to  detect  any  y-rays  from  a  CSO. 

The  LAT  radio  galaxies  and  unclassified  objects  are  a  puzzle. 
They  do  not  appear  to  differ  significantly  from  their  non-LAT 
counterparts  in  any  appreciable  way.  Their  opening  angles 
appear  to  be  very  similar  to  the  non-LAT  objects  and  they  are 
not  particularly  bright  in  the  radio.  The  LAT  radio  galaxies/ 
unclassified  objects  may  have  higher  core  temperatures  than 
the  non-LAT  objects.  The  one  difference  that  stands  out  is  that 
we  detected  polarization  in  the  cores  of  67%  of  the  LAT  radio 
galaxies/unclassified  objects  compared  to  only  20%  for  their 
non-LAT  counterparts.  Unfortunately,  our  current  sample  is  too 
small  to  make  any  definite  claims  about  these  objects.  There  is 
also  controversy  about  the  optical  types  of  several  of  the  objects 
we  label  as  radio  galaxies. 

6.4.  Constraints  on  the  y-ray  Emitting  Region 

The  major  differences  between  LAT  and  non-LAT  sources  are 
confined  to  the  cores  of  the  sources.  The  LAT  BL  Lac  objects 
have  stronger  jets,  suggesting  higher  ejection  velocities  from  the 
cores.  The  LAT  FSRQs  have  higher  core  temperatures  and  more 
frequently  exhibit  polarized  cores.  The  LAT  radio  galaxies/ 
unclassified  objects  are  also  more  polarized  than  their  non-LAT 
counterparts.  We  suspect  that  the  y-rays  are  being  produced 
very  near  to  the  cores  of  these  objects.  To  put  a  limit  on  the 
size  of  the  y-ray  emitting  region,  we  measured  the  sizes  of 
the  cores  for  several  objects.  Using  all  of  the  objects  for  which 
we  have  measured  redshifts,  the  median  upper  bound  on  the  core 
size  is  4.54  pc  with  a  standard  deviation  of  10.98  pc.  However, 
a  better  estimate  for  the  upper  limit  on  the  size  of  the  emitting 
region  should  come  from  measurements  of  the  objects  nearest 
to  us.  Limiting  ourselves  to  objects  with  z  0.1,  we  have  four 
BL  Lac  objects  and  three  radio  galaxies  (a  LINER,  a  Seyfert  1, 
and  an  object  that  may  actually  be  an  FSRQ).  Using  this  subset 
of  sources,  we  find  a  median  upper  bound  on  the  core  size 
of  0.58  pc  with  a  standard  deviation  of  0.34  pc.  Therefore,  the 
y-ray  emission  region  must  be  within  about  0.9  pc  of  the  central 
engine. 

7.  CONCLUSIONS 

We  compared  the  distributions  of  core  polarization,  core 
brightness  temperature,  jet  opening  angles,  core  fractions,  jet 
bending,  jet  length,  and  8.5  GHz  radio  flux  densities  for  LAT 
and  non-LAT  sources.  The  only  significant  difference  between 
the  LAT  and  non-LAT  BL  Lac  objects  is  that  the  y-ray  loud 
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BL  Lac  objects  appear  to  have  longer  jets.  Also,  we  find  no 
correlation  between  BL  Lac  radio  flux  density  and  y-ray  flux. 
We  suspect  that  all  BL  Lac  objects  are  producing  y- rays,  but 
some  are  below  LAT’s  threshold  due  to  jet  material  speed,  or 
possibly  variability  in  the  y- ray  luminosity.  In  contrast,  y- ray 
bright  FSRQs  are  considerably  different  from  the  y- ray  quiet 
population.  The  LAT  FSRQs  appear  to  be  extreme  sources  with 
higher  core  brightness  temperatures,  larger  opening  angles,  and 
greater  core  polarization  than  their  non-LAT  counterparts.  There 
is  also  a  significant  correlation  between  the  FSRQ  radio  flux 
density  and  y- ray  flux,  especially  for  brighter  sources.  The 
y- ray  loud  FSRQs  can  be  explained  by  Doppler  boosting,  but 
the  orientation  and/or  jet  material  speeds  must  be  significantly 
different  than  the  y-ray  quiet  FSRQs.  One  possible  way  to 
explain  this  difference  is  that  the  y-ray  loud  FSRQs  may  have 
greater  angular  momentum  in  their  central  supermassive  black 
holes,  leading  to  higher  jet  velocities  and  greater  polarization. 
Due  to  the  small  number  of  radio  galaxies  in  our  sample, 
combined  with  the  controversy  over  some  of  their  optical  types, 
there  is  little  we  can  say  about  y-ray  bright  radio  galaxies. 
However,  we  do  find  two  intriguing  clues  about  them.  First, 
there  are  no  y-ray  loud  CSOs  in  our  sample.  Second,  the  radio 
galaxies  that  we  do  detect  in  the  y- rays  tend  to  have  strong 
polarizations  in  their  cores. 

Because  the  differences  we  find  between  the  y-ray  loud  and 
y-ray  quiet  objects  are  related  to  the  core,  we  suspect  that  the 
y-ray  radiation  originates  within  the  core  (i.e.,  at  the  base  of 
the  jet).  By  calculating  the  core  sizes  of  the  LAT  objects  with 
measured  redshifts,  we  put  an  upper  limit  of  15  pc  on  the  size  of 
this  emitting  region.  A  better  limit  of  0.9  pc  results  from  using 
the  core  sizes  of  the  lowest  redshift  objects  in  our  LAT  sample. 

As  more  data  are  collected  more  will  become  clear  about 
y-ray  loud  blazars.  We,  and  several  other  groups,  are  currently 
conducting  further  observations.  We  will  soon  have  more 
information  about  the  velocities  of  the  jet  material  and  more 
associations  of  AGNs  with  y-ray  sources. 
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